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We investigate the sedimentation of knotted polymers by means of the stochastic rotation dynamics, a molec- 
ular dynamics algorithm which takes hydrodynamics fully into account. We show that the sedimentation co- 
efficient s, related to the terminal velocity of the knotted polymers, increases linearly with the average cross- 
ing number n c of the corresponding ideal knot. To the best of our knowledge, this provides the first direct 
computational confirmation of this relation, postulated on the basis of experiments in Ref. [1], for the case of 
sedimentation. Such a relation was previously shown to hold with simulations for knot electrophoresis. We 
also show that there is an accurate linear dependence of s on the inverse of the radius of gyration Rg 1 , more 
specifically with the inverse of the R g component that is perpendicular to the direction along which the polymer 
sediments. Intriguingly, the linear dependence of s on n c notably deteriorates in the presence of walls, even 
when the polymer sediments in a slab which is over an order of magnitude wider than the typical size of the 
knots we study. In contrast, R~ x remains to a good precision linearly dependent on n c . Therefore, R~ x is a 
good measure of a knot's complexity. 



PACS numbers: 87.15.A-,82.35.Lr,82.37.-j 

I. INTRODUCTION 

The topology of knots has been studied for almost two cen- 
turies. The earliest work dates back to the 19th century [2, 3], 
when P. G. Tait first proposed the current classification of 
knotted loops based on the crossings observed in their two- 
dimensional projection. Since then, knot theory has been ap- 
plied to a wide number of areas in physics and beyond: of 
particular relevance to our work are the applications to bio- 
physics, as several biopolymers are observed to form knots. 
A particularly important case is that of DNA [4-10]. DNA 
has been long known to routinely form simple knots in solu- 
tions [11], and it has been argued that specific enzymes must 
be at work within a cell to avoid excessive DNA knotting [12], 
which could be detrimental for the correct genome function. 
More recently, knots have been discovered also in the native 
states of proteins, where they are however quite rare [13]. 
Knots are further present in synthetic polymers. 

DNA knots are even more common within bacteriophage 
heads [8], where they form due to the tight geometric confine- 
ment the genome is subjected to. Such bacteriophages are the 
"viruses of bacteria", and they work by releasing their DNA 
into the cytosol of their host. This release, or DNA ejection, is 
essentially entropically driven, as the DNA of phages is con- 
fined to essentially crystalline density, which is highly costly 
entropically. Also wildtype phages when burst open reveal a 
knotted double- stranded DNA. Modelling has shown that the 
spectrum of knots observed in the bacteriophage P4 [8] can 
be understood on the basis of simple polymer physics models 
that view DNA as a semiflexible and self-avoiding polymer, 
provided that an aligning interaction between contacting DNA 
segments is included [7]. 

As the knots are so frequent and important in the biophysics 
of DNA, there is clearly a demand for their easy identification. 
But how can one tell whether a DNA molecule, whose (hy- 
drated) thickness is about 2.5 nm, is knotted, and which knot 
it forms? In a seminal paper [4], Stasiak and co-workers pre- 



sented experimental results on DNA gel electrophoresis and 
showed that the measured electrophoretic mobility v of the 
DNA knots used in the experiments increases linearly with the 
average crossing numbers n c of the ideal forms (defined be- 
low) of these knots [4, 5]. This dependence makes the identi- 
fication of DNA knot topology much faster than e.g. by using 
electron microscopy methods; it also allows the determina- 
tion of knotting probability in an ensemble of DNA molecules 
without the need to examine each of these singularly. 

How can we compute n c in practice? To do so, we first need 
to deform a knot into its ideal form, which is defined as the 
one with the highest ratio of volume to surface area [14, 15]. 
In other words, this is a knot that is formed with as short a 
ropelength (polymer) as possible. The average crossing num- 
ber, n c , is defined as the average number of crossings of all 
possible two dimensional projections of the knot in its ideal 
form [14]. Clearly, the larger n c becomes, the more complex 
is the corresponding knot. Knots are typically indexed accord- 
ing to their minimal number of crossing, n m i n , which is the 
minimum number of crossings over all projections - therefore 
n c > fhnin.Le. the number of crossings that cannot be opened 
without breaking the knot contour. 

The reason why the electrophoretic mobility of knots in- 
creases with average crossing number, hence their complexity, 
can be understood with a simple physical argument. If a DNA 
knot is subject to a force /, for instance due to an electric field, 
then the (terminal) velocity it will reach can be estimated by 
the formula / = jv, where 7 is an effective friction - this 
is because the dynamics of DNA is highly overdamped in so- 
lution. A simple guess for 7 is then given by Stokes law as 
67rr]Rg where 77 is the solvent viscosity, and R g is the radius 
of gyration of the polymer, which measures its typical "size". 
If the knots have the same contour length, it is then intuitive 
that the more complicated ones will be more globular, hence 
smaller in shape: their terminal velocity, hence mobility, will 
then be larger. 

This line of reasoning is essentially the one in [1, 4], where 
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the authors used a method to calculate the expected sedimen- 
tation coefficient s of DNA molecules with a given topol- 
ogy [1], and found that s increases linearly with n c . How- 
ever, the original computation that employed the Kirkwood- 
Riseman approximation neglected the effect of flow on poly- 
mer conformation, and was essentially an estimate only valid 
at infinitesimally small forcing. To quantitatively establish the 
claim on the linear relationship between v or s and n c one 
needs to simulate the sedimentation of knotted polymers us- 
ing a computational method that studies the molecular dynam- 
ics of the polymer in the presence of hydrodynamics and of a 
gravity field. This is what we set out to do in the present paper. 
As our work comprises the first numerical simulations of knot 
sedimentation including full hydrodynamic interactions, our 
results and framework may be used in the future to enhance 
accuracy in knot determination by sedimentation experiments. 
We also highlight the importance of boundaries, which signif- 
icantly affect the sedimentation coefficients we record. Our 
results may be seen as complementary to the dynamic Monte- 
Carlo simulations in Ref. [16, 17] which established the lin- 
earity of the electrophoretic motility with n c via direct sim- 
ulations of the dynamics of a knot in a gel, under an electric 
field. 

Our paper is organized as follows. In Sec. II we explain our 
computational model for polymer sedimentation: the model 
for the DNA, the method for simulating the dynamics and the 
used simulation geometries. In Sec. Ill we present the results 
and conclude in Sec. IV. 



II. THE COMPUTATIONAL MODEL 

A. The polymer model 

The polymer is modeled as a circular chain of pointlike 
beads with mass ra^. The adjacent pairs are connected by 
means of the finitely extensible nonlinear elastic (FENE) po- 
tential 

Up = -yi? 2 ln(l - (^) 2 ) , r<R, (1) 

where r is the length of the bond and R = 1.5<j is the max- 
imum bond length. All bead pairs interact through a shifted 
truncated Lennard- Jones (LJ) potential, 




(2) 

where is the distance between beads i and j. The po- 
tential parameters are chosen as a = 1.0, e = 1.0, and 
K — 30/a 2 . This truncated repulsive Lennard- Jones poten- 
tial models a good solvent. 

We also examine the effect of a simple bending potential 
model where every consecutive segment pair feels a potential 
of the form 

^bend = ~K (r i+ i ~ Ti) • [fi ~ fi-i) . (3) 



Here, fi is the location of bead i. We run the simulations with 
k = 0, corresponding to a flexible polymer, and with k = 2, 
corresponding to a semiflexible polymer, or worm-like chain 
(WLC). 



B. Polymer and solvent dynamics 

In the modeled sedimentation process the polymer is im- 
mersed in a solvent and is driven by a constant (gravita- 
tional) force. This makes sedimentation an inherently non- 
equilibrium process. [For such processes hydrodynamics typ- 
ically plays a significant role.] Hydrodynamic interactions are 
taken into account by modeling the solvent using stochastic 
rotation dynamics (SRD), a computationally efficient Navier- 
Stokes integrator [18]. We use a hybrid version of the al- 
gorithm where the polymer follows Newton's dynamics and 
SRD is applied to both the solvent and the polymer. The poly- 
mer molecular dynamics is implemented using the standard 
velocity Verlet algorithm [19]. The above described model 
has been succesfully used in modeling elasticity and hydrody- 
namics of linear DNA molecules [20]. 

The solvent consists of pointlike particles, each of mass m s . 
One SRD step comprises two smaller steps: the streaming step 
and the collision step. In the streaming step the positions of 
solvent particles are updated as 

fi(t + At) = fi(t) + Vi(t)At, (4) 

where fi (t) and fi (t) are the location and the velocity of the 
particle i at time t, respectively, and At is the SRD timestep. 
The whole simulation space is divided into a grid of cubic 
cells of equal size. In each cell the velocities of particles are 
updated in the collision step as 

Vi(t + At) = v cm (t) + R [vi(t) - v cm (t)] , (5) 

where v cm (t) is the center-of-mass velocity of the cell. Oper- 
ator R is a rotation matrix whose rotation angle a is fixed but 
the rotation axis is chosen randomly in each timestep for each 
cell. The method conserves energy and momentum in each 
cell. In order for the method to maintain Galilean invariance 
the grid is shifted randomly at each step [21]. 

The polymer is coupled to the solvent in the collision step, 
Eq. (5). Here, the polymer beads are treated like the solvent 
particles. The resulting polymer beads' velocities are then 
used in the following velocity Verlet steps performing molec- 
ular dynamics (MD). MD and SRD take turns so that after 
every 500 MD steps of timestep St = 0.002 is performed a 
single SRD step of timestep At = 1. The other parameters 
were chosen as follows: the edges of the cubic cells are of 
unit length, the average density of solvent particles is 5 parti- 
cles per unit cell, and the rotation angle a = 37r/4. The mass 
of the solvent particles is m s = 4 and the mass of the polymer 
beads is ra& = 16. 
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C. The simulation geometry and polymer configurations 

The simulation space is a parallelepiped with periodic 
boundary conditions along the y and z directions for both 
the polymer and the solvent particles. Polymers are driven 
by a constant force in the negative z direction. To model dif- 
ferent experimental setups where sedimentation takes place 
in a large container and in a channel we run the simulations 
with either periodic boundary conditions also along the x di- 
rection, or in the presence of two no- slip walls parallel to 
the yz plane and normal to the x direction. When the sol- 
vent particle hits the no-slip wall its momentum is reversed, 
so it is bounced back to the direction of incidence [22]. In 
addition to this there is a repulsive Lennard- Jones potential 
similar to Eq. (2) with a cutoff at distance 0.1 (in simulation 
units) from the wall. The dimensions of the parallelepiped are 
L x x L y x L z = 50 x 50 x 300 when the periodic boundary 
conditions are applied and L x x L y x L z = 10 x 50 x 300 
when the no- slip walls are present. 

In each simulation we use a circular polymer with the num- 
ber of beads 7V& = 215. However, initial configurations dif- 
fer since we examine polymers with different knot topologies. 
The different knot topologies are created with the Knotplot 
program. Individual realizations of the real conformations 
corresponding to a single knot topology were created by ther- 
malizing the knots using different stochastic forces due ther- 
mal fluctations. Snapshots of the initial configurations and 
of typical conformations (during sedimentation) are shown in 
Fig. 1, for three different knot types. 

On each polymer bead we apply a constant force / = 

a = 0.1 in the negative z direction, where a is the acceler- 
ation. Once the polymer has reached its terminal velocity v t 
we calculate the time average of the sedimentation coefficient 



in the negative z direction. Applying the constant force to the 
polymer induces momentum in the polymer which transfers 
to the solvent via the SRD collision step, Eq. (5). 

When periodic boundary conditions are applied in the z- 
direction the unphysical increase in the momentum of the fluid 
that in effect rotates in the z-direction has to be judiciously 
taken into account and compensated for. This we do by re- 
moving the momentum I/N s from every solvent particle at 
each streaming step. I = fN^At is the impulse applied to the 
polymer due to the constant force between SRD-steps. The 
same impulse is then divided equally to each of the N s sol- 
vent particles. The procedure preserves the local momenta in 
the fluid induced by the sedimentating polymer but removes 
the unphysical net momentum building up due to the periodic 
boundaries. 

The no- slip walls instead naturally dissipate the induced en- 
ergy. Hence, there is no need for momentum correction for the 
sedimentation simulations within the channel. For both cases 
the total momenta of the solvent in the x— and directions 
are kept at zero. The solvent is kept at a constant tempera- 
ture by scaling the momenta so that the equipartition theorem 
always holds. 



(A) (B) (C) 




FIG. 1. (Color online) Initial configurations (top row) and typi- 
cal configurations during sedimentation runs (bottom row) for three 
knots: trefoil (3i, (A)), Stevedore's knot (6i, (B)) and 10i knot (C). 
[These are three examples of twist knots, with the trefoil being a 
torus knot as well.] 



III. RESULTS 

The sedimentation coefficient for each polymer knot topol- 
ogy was obtained by averaging over 200 sedimentation simu- 
lations. The standard error of the mean was used as the esti- 
mate of error. In what follows, we denote an affine dependen- 
dence of a variable y on a variable x, i.e. y = Ax + B, where 
A and B are constants, as y ~ x. Also, by 'linear depen- 
dence' we mean affine dependence. [Stasiak and coworkers 
use linear dependence in this same manner.] 

First we investigate the sedimentation of flexible knots 
(k = 0) in the bulk (i.e. by using periodic boundary conditions 
rather than confining walls). Fig. 1 shows initial and typical 
configurations during sedimentation runs for three knots. In 
Fig. 2(a) the sedimentation coefficients s for polymer knots 
are plotted against the average crossing numbers n c for their 
respective ideal forms. These crossing numbers are obtained 
from [15]. A very good linear relationship between s and n c 
is obtained for the common DNA knot topologies that are also 
reported in the original papers of Stasiak and coworkers. For 
the more complicated, higher order, knots, the dependence of 
s on n c deviates somewhat from linearity. Modulo this, our di- 
rect molecular dynamic simulations (including hydrodynamic 
interactions fully) confirm well that there should be a linear 
increase of s with n c , as proposed by Stasiak and coworkers, 
at least for relatively simple knots. 

Stasiak and coworkers explained the linear relationship be- 
tween n c and the knot migration speed as due to the fact that 
n c is directly proportional to the compactness of the knots. As 
a measure of molecular compactness they used the mean of in- 
verse distances in a molecule [23], [R -1 ] = ^fj^/j)^' 1 )' 
where (r^ 1 ) is the mean reciprocal separation of segments i 
and j [4]. The quantity [R~ r ] indeed appears in the approxi- 
mate hydrodynamic theory used to compute sedimentation co- 
efficients from polymer conformations [1] - at the same time 
though this is hard to calculate in practice (see e.g. [24]). 
Hence, it is not easy to show that n c increases linearly with 
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FIG. 2. (Color online) Periodic boundary conditions, (a) Sedimen- 
tation coefficient s and (b) inverse of the radius of gyration Rg 1 as 
a function of knot average crossing numbers n c for respective ideal 
knot conformations. Knots from left to right in all figures as function 
of n c are as follows: 3i, 4i, 5i, 52, 61, 62, 63, 7i, 73, 72, 7s, 74, 76, 
7 7 , 81, and9i. 
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FIG. 3. (Color online) Periodic boundary conditions. Sedimentation 
coefficient s for knots with node numbers (a) rid — Q and (b) rid — 7 
as functions of knot average crossing number n c . Within the knot 
family rid = 7 the linear relation s ~ n c does not hold. 



[R~ x } as would be required for the linear relationship between 
n c and v to be proved rigorously. 

The inverse radius of gyration, i?~ 1 is also, by definition, 
a good measure of molecular compactness, and it can be 
readily and efficiently computed in simulations, as R g = 

(1/7V) Xli=i(( r i — Rem) 2 ), where R cm is the center of mass. 
In Fig. 2(b) n c is shown as a function of the inverse radii of 
gyration measured for the sedimenting polymer knots. We ob- 
tain a good linear relation Rg 1 ~ n c for the knots investigated 
by Stasiak and coworkers. Interestingly, the slight deviations 
from the linear dependence s ~ n c show up in almost ex- 
actly the same way in Rg 1 ~ n c . [That is, also the relation 
R~ x ~ n c tends to break down for highly complex knots such 
as 81 and 9i - just as the s ~ n c relation did.] 

By a careful look at the s vs n c curves, one may note that 
the sedimentation coefficients for the knots within knot fam- 
ilies with 6 and 7 minimal crossing deviate from the linear 
dependence s ~ n c . The sedimentation coefficients s ob- 
tained for all knots with n m i n = 6 and n m i n = 7 are shown 
in Fig. 3. For the knots 61, 62, and 63, s has a close to lin- 
ear dependence on n c - however, the slope differs from the 
one characterising the general linear dependence on n c for all 
the knots. Within the knot family with n m i n = 7, s violates 
linearity more significantly. For more complex knots (such as 
81 and 9i), deviations from linearity increase, as previously 
mentioned. All this is in good accord with the experimental 
observation that the linear relationship between the speed of 
electrophoretic migration speed and n c breaks down for com- 



FIG. 4. (Color online) Sedimentation coefficient s as a function of 
(a) inverse radius of gyration Rg 1 and (b) inverse of perpendicular 
component of radius of gyration R~^ . In (a) the dashed (red) line is 
the estimate for s, Eq. (7), based on the Stokes' formula. Here, only 
the knots with subscript 1 are depicted Also the unknot 0i and the 
knot 10i are included. Periodic boundary conditions were used to 
generate all these data. 

plex knots [25]. 

Fig. 4 (a) shows s as a function of R g for the first knot type 
(subscript 1) of each knot family n m i n G {0,3,4,... ,10}. A 
very accurate relation s = ARg 1 + B, where A = 2.16 and 
B = 1.1, is obtained. In other words, the linear proportional- 
ity between s and Rg 1 (at least for sedimentation runs in the 
bulk) appears to hold more accurately than the one between s 
and n c . [That is because the dependencies of s and Rg 1 on 
n c show similar deviations from linearity, as observed above.] 

The dependence s ~ R~ x is in agreement with the simple 
physical argument given in the introduction; it has also been 
derived for DNA supercoils [26], based on the Kirkwood- 
Riseman expression valid for linear Gaussian chains of radius 
of gyration R g [27]. Using the Kirkwood-Riseman expres- 
sion entails assuming that the sedimentation is non-draining, 
i.e. that the motion of the solvent particles in the region of 
the polymer is largely suppressed (this is essentially the same 
approximation used in computing sedimentation coefficients 
in [1]). Consequently, backflow couplings between solvent 
flow and polymer conformations are completely neglected in 
that approach (unlike in our SRD simulations). 

It is instructive to compare the values of s measured from 
our simulations with the result expected from the Stokes' for- 
mula for friction (see Introduction), which would be 



This estimate is depicted by the dashed line in Fig. 4 (a). Here, 
M is the total mass of the polymer and r] = 4.78 is the vis- 
cosity of the fluid estimated for SRD using formulas derived 
in Ref. [28]. This rather crude estimate, which views the knot 
as a sedimenting solid ball of radius R g and mass M, appears 
to describe our data reasonably well, considering the level of 
the approximation. The Stokes estimate is most precise for 
the knots 5i, 61, and 7i. 

Even though these simplified theories (Stokes drag and 
Kirkwood-Riseman non-draining approximations) work rea- 
sonably well, the effect of the velocity field on the polymer 
(backflow) cannot be wholly disregarded (nor can the near- 
field flow which is not captured by such approximations). The 
most relevant effect in our simulations which is not captured 
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FIG. 5. (Color online) (a) Different, appropriately normalized, com- 
ponents of radius of gyration with periodic boundary conditions. 
Spheres: R g — R g /V3, squares: R g ± — R g ±/V%, triangles: R g \\. 
(b) Different cartesian components of radius of gyration with slit 
walls. Sphere: R gx , square: R gy , triangle: R g \\ = R gz . 
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FIG. 7. (Color online) Results for sedimentation within a slab, with 
no-slip walls parallel to the direction of sedimentation. Sedimen- 
tation coefficient s as a function of (a) the inverse of the radius of 
gyration R g and (b) the inverse the perpendicular component of the 
radius of gyration R g ± 
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FIG. 6. (Color online) Slit walls, (a) Sedimentation coefficient s as 
a function of average crossing number n c . (b) Inverse of average 
radius of gyration Rg 1 as function of knot average crossing number 
n c . 



in the more simplified treatments is that the components of 
R g aligned, R g u, and perpendicular, R g ±, to the direction of 
gravitational force differ appreciably, see Fig. 5 (a), at least 
for the values of a used here. The dominant contribution to 
the linear dependence s ~ R- 1 comes from s ~ Rg~±> see 
Fig. 4 (b). Both and R~]_ depend linearly on n c to a 
good precision. However, the linear dependence s ~ 
(not shown) breaks down for knots of higher complexity, 9i 
and 10 1. This in keeping with the friction being determined 
dominantly by the dimension R g ±. 

In order to consider the effect of boundaries, which are not 
necessarily negligible in the lab (or could be used to control 
the drag exerted by the fluid on the polymer), we also simu- 
lated sedimentation of knotted polymers between two no- slip 
walls. Due to the interaction of the fluid particles with the 
walls the friction felt by the polymer is enhanced, and as a 
result the polymer conformation is more elongated in the di- 
rection parallel to the walls. Hence, R g \\ deviates from R g ± 
more clearly in this case with boundaries. Interestingly, now 
the linear relation between s and n c is not as precise as in 
the bulk (periodic boundary conditions, without walls), see 
Fig. 6 (a). This is explained by Rg 1 deviating slightly more 
from its linear dependence on n c , see Fig. 6 (b). Accordingly, 
the relation s ~ Rg 1 in the channel, Fig. 7 (a) is far from the 
precise relation obtained with periodic boundary conditions. 
From Fig. 5 (b) it is clear that the component of the radius of 
gyration that is perpendicular to the walls, R gx , is squeezed 
and has no dependence on n c . The width of the channel is 8 
and large compared to R g « 0.7 for the trefoil in free solvent. 



FIG. 8. (Color online) Sedimentation coefficient s (circle) and radius 
of gyration R g (triangle) of the knot 6i as a function of the bending 
potential parameter k. Periodic boundary conditions. 

Hence, the precision in identifying knot topologies through s 
deteriorates quite dramatically under confinement. However, 
the linear dependence of s on R~y , the component measured 
in the perpendicular direction on which periodic boundaries 
are applied, is preserved to a good precision, see Fig. 7 (b). 

Adding the bending potential Eq. (3), for k ^ 0, R g in- 
creases with increasing k. The sedimentation coefficient s 
diminishes accordingly with increasing k, see Fig. 8. Sedi- 
mentation runs for worm like chains (WLC) (k ^ 0) leads 
to the same qualitative trends for s as a function of n c as ob- 
served for flexible chains (k = 0). However, the deviations 
from the affine relationships observed for the flexible chains 
appear to quantitatively increase with increasing k. Also, the 
linear dependence of s on Rg 1 for WLC deviates from the 
estimate given by Eq. (7) more than this dependence for flex- 
ible chains. This is in accord with what was stated in con- 
nection with Fig. 4 (a): a WLC sediments in a sparser (less 
compact) conformation than a flexible chain, hence, Stokes' 
formula will describe its friction less well. 



IV. CONCLUSION 

In conclusion we have investigated sedimentation of knot- 
ted polymers using stochastic rotational dynamics, a compu- 
tational model where hydrodynamics is taken fully into ac- 
count, for the first time. Our motivation was to directly test 
the dependence of the sedimentation coefficient s on the ideal 
average crossing number n c , without making any assumptions 
about the segment distribution or hydrodynamic interactions. 
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The linear dependence of s on n c was seen to hold with fair 
to good accuracy for the knot topologies for which Stasiak 
and co-workers predicted this dependence. However, for some 
knots of higher complexity (such as 7-crossing number knots, 
81 and 9i), deviations from this dependence increase. This 
is in line with the experimental and theoretical results ob- 
tained for gel electrophoresis of knots [16, 17, 25]. Our direct 
sedimentation simulations further justifies the argumentation 
commonly used in the literature to explain the observed lin- 
earity of s on n c in gel electrophoresis experiments. 

Our simulations show that the inverse radius of gyration of 
the polymer, Rj 1 , is also proportional to the average crossing 
number n c , and therefore provides another good measure of a 
knot's complexity. Interestingly, the deviations from linearity 
in the R~ x vs n c curve pretty much mirror those in the s vs 
n c curve. As a result, we observe that the linear relationship 
between s and R~ x is far more precise than the one between 
s and n c . As the solvent flow affects the knot conformations, 
and renders them more anisotropic, we monitored separately 
the dependence on the component of the radius of gyration 
parallel and perpendicular to the sedimentation direction: we 
found that it is the inverse of the latter which is more accu- 
rately linearly proportional to s. 



Our fluctuating hydrodynamics simulations finally suggest 
that confinement, or the presence of boundaries, may signifi- 
cantly affect the results, and lead to a deterioration of the lin- 
earity between s and n c . Even when the distance between the 
walls is an order of magnitude larger than the radius of gyra- 
tion of the polymer (in the bulk), these effects are quite large. 
The deviations from linearity in the s vs R~ x curves are due 
to R g component perpendicular to the walls losing all depen- 
dence on n c . 

We hope that our simulations will spur further results on 
knot sedimentation, aimed at verifying the deviations from 
linearity in the s vs n c curve which can be observed due to 
confinement, or at high knot complexity. Eventually, find- 
ing a more precise relation for complex knots will lead to an 
improvement in the knot identification techniques when the 
average crossing number is high. 
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